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Page L, line 3' Replace the quantity - 1 by \1 m^ - 1- 

Page 11, equation (17): In the first value for |i- for plan form D, 

T ^ 1 

the term — within the final psu:*entheBeB should be — . 

2b 2b 

Page 21: Correct equations (B3) as follows: 

In line 2, the expression - Ep-| should be a^Epi . 

In lines k, 6, 7, 13, I5, and I8, all primes appearing on the right- 
hand side of these equations should be replaced by the superscript 1. 

In line li<-, the sixth term of the equation for -2J1M5' should be 
5 jHi 2® instead of 5^Hip. 
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Page 23: In lines 2, 3s and 4, the equations for Bg, By, and B3 
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By Herbert C. Nelson 


SUMMARY 


Expressions for the velocity potential and the lift and moment are 
derived for thin, oscillating, "arrowhead" type wings with supersonic 
edges moving at a constant supersonic speed. The triangular wing and 
any one of the series of wings obtained from it by cutting its rearward 
part so that the flow normal to the resulting trailing edges remains 
supersonic are Included in the term arrowhead type. Explicit results 
for the lift and moment are given for the wings undergoing harmonic 
pitching and vertical oscillations. 

Closed expressions for the velocity potential, section force and 
moment coefficients for any arrowhead wing, and total force and moment 
coefficients for only the triangular or delta wing are developed explic- 
itly to the third power of the frequency of oscillation. These expres- 
sions may be sufficient for most practical applications. 

The wings considered are found to exhibit the possibility of nega- 
tive damping in torsion for certain ranges of Mach number and axis-of- 
rotatlon position. A figure showing the ranges of these parameters for 
triangular wings is given herein. 


INTRODUCTION 


In high- speed-aircraft design, a knowledge of the air forces and 
moments that act on various types of oscillating-wing plan forms is 
often required. Such knowledge' is useful in .considering a number of 
instability problems, among which are wing flutter and low-frequency 
instability of the aircraft involving control- surface deflections. The 
main source of theoretical information has been the solution of the 
linearized differential equation for compressible flow. 

The problem of finding the air forces on an oscillating or steady 
thin wing in supersonic flow is us\aally formulated in terms of the dis- 
turbance velocity potential for the wing. The disturbance velocity 
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potential is constructed ty superposing solutions of the linearized 
differential equation, corresponding to sources, doublets, or higher- 
order singularities, in such a way as to satisfy the boundary conditions 
at the Wing surface and at infinity. The linearized differential equa- 
tion and its boundary conditions constitute the boundary- value problem 
for the velocity potential for the wing. 

The present paper is. concerned with the boundary- value problem for 
the velocity potential for an apex-forward, thin, flat, oscillating, 
triangular or related wing with supersonic edges moving at- a constant 
supersonic, speed. According to reference 1, this problem may be clas- 
sified as "purely supersonic," since the upper and lower surfaces of the 
wing can be regarded as acting independently of— one another, and is sat- 
isfied by a surface distribution of souircea with local strength propor- 
tional to the local prescribed normal velocity at the wing surface. 
Although reference 1 does hot explicity trelt the oscillating triangular 
wing^ its general solution must be regarded as giving the velocity 
potential for this wing in integral form, ^he purpose of the present 
paper is then to obtain an integrated . form for the. .velocity potential, 
from the integral-form solution of reference. 1, for the particular case 
of the oscillating triangular or related wing. 

In a strict sense this integrated form Applies to a semi -infinite 
triangular plan form. But-^ since the wing imke is o:g no concern in the 
present development, the Velocity potential_may be considered to apply 
to any one of the series of wings obtained from a triangular wing with 
supersonic edges by cutting the rearward part of this wing in such a 
manner that the flow normal to the resulting trailing edges remains 
supersonic. All these wings, including the . triangular wing, are hence- 
forth referred to in general as arrowhead wings. 

Thete ate several other papers closely associated with the subject 
being considered. In references 2 and 3 the treatment of oscillating, 
finite, s^pt wings involved an' Integral form for the velocity potential 
quite similar to that for the oscillating arrowhead wing. In refer- 
ences 4 and more directly, eicpressions for the total lift and moment 
were given for the particular oscillating arrowhead wing with super- 
sonic edges known as the "wide'' delta wing. 

The present paper effectively includes "the total forces and moments 
for the wide delta given in references 4 and 5* In addition the' present 
paper gives the forces and moments on any streamwlse wing section for 
not only the wide delta bilt also for the mote general arrowhead wing. 

The section rather than the total forces and moments are desirable, for 
example, in a strip flutter analysis that includes the wing-flutter mode 
shape. 

The velocity potential for the arrowhead wing does not appear to 
be obtainable in terms of elementary functions (see references 1 and 4 ). 
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The solution of this type of prohlem, therefore, is approximated hy 
various kinds of expansion. Thus in reference the velocity potential 
was expanded in a series of Bessel functions; whereas in reference 6, 
in the treatment of the rectangular wing, an expansion in powers of the 
frequency of oscillation was employed. The frequency-expansion method 
is utilized in this paper and only the first few terms of the velocity 
potential and subsequent forces and moments are obtained. The first 
few terms are considered adeqiiate for a large class of practical appli- 
cations. In the derivation of the velocity potential and the forces 
and moments, the wings are assumed to be undergoir^ harmonic oscilla- 
tions in pitch and vertical translation. 


SYMBOLS 


(fi disturbance-velocity potential 

x,y, z . rectangular Cartesian coordinates 

U velocity of main stream 

a velocity of sound in main stream 

t time • 


w(x,y,t) normal velocity at surface of wing at point (x,y) 

Zjj^ function defining vertical displacement of point (x,y) 

of wing 

xq abscissa of axis of rotation of wing section as shown in 

figure 1 

p air density in main stream 

h vertical displacement of axis of rotation 

h^ maximum amplitude of vertical displacement of axis of 

rotation^, positive downward 

a angle of attack 

ttQ maximum amplitude of angular displacement about axis of 

rotation, positive leading edge up 

• • 

h,a time derivatives of 


h and a, respectively 
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CD 


circular frequency of oscillation 


M free-stream Mach number (u/a) 

3 = - 1 




CD 



®'n^^m 


rectangular Cartesian coordinates used to represent 
location of sources in xy-plane 


functions of cd, x, and M defined in equation (ll) 


X 


slope of leading edge of wing as shown in figure 1 


■5 = lim 
^X—^1 


u,v characteristic coordinates defined in appendix A 


. -LjlM 
' 1 + 3X 

2 2 

(J = 3 4^ - 1 

b semichord of midspan wing section 

|i.,V^lJ.Q nondimensional coordinates 

k reduced frequency (bcD/u) 

Ap local- surface pressirre difference 



P 


Li, Mi 


M, 


a 


section force per unit span on wing strip parallel to main 
stream, positive downward 

section moment per unit span on wing strip parallel to main 
stream, taken about axis of rotation Xq, positive leading 
edge up 

components of section force and moment coefficients, 
respectively, defined in equation (B1 )j i = 1, 2,3,4 

total force on delta wing, positive downward 

total moment'on delta wing about axis of rotation x = Xq, 
positive leading edge up 
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components of total force and moment coefficients, 
respectively, defined in equation (bU); i = 1,2, 3, 4 


I? ■R' 
•“mn'> mn-> 

G ,H 
mn mn 

“i^PiT] 

7i,8iJ 


functions of 


functions of 


3 , and V defined in equation (B 2 ) 
M, and 0 defined in equation (B3) 


Z 


s 


distance to point at which extended trailing edge of plan 
form D, figure 3^ intersects x-axis 

semi span of wing 


FORMULATION OF PROBLEM 


Consider uniform supersonic flow past a thin, flat, arrowhead type 
wing as shown in figure 1 (with its leading edges outside the Mach cone 
generated by the nose of the wing) . The wing is referred to a coordinate 
system fixed in space and is assumed to he creating small disturbances 
in the main stream flowing past. Then, if the undisturbed stream 
velocity is in the direction of the positive x-axis and is of magni- 
tude U, the equation satisfied by the distvirbance velocity potential 
for the wing is 

^(it (1) 

The wing is assumed to be performing oscillations of small ampli- 
tude about its own undisturbed position, the plane Z = 0 . The boundary 
condition to be satisfied at the wing surface is then 

where is the vertical displacement of a point (x,y) of the wing. 

Note that this boundary condition is evaluated, in accordance with 
small- disturbance linearized theory, at the xy-projectlon of the wing 
rather than at the wing Itself. The additional boundary conditions, 
that only the wing can support a pressure difference and that the sources 
of disturbance must not be felt ahead of their respective Mach cones, 
are automatically satisfied by the type of source synthesis to be 
employed in the solutions. 
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After the boundary-value problem is solved for the velocity 
potential the pressure on either wing sxjrface may be found by means 
of a linearized form of Bernoulli's equation 

p - D 

where p is the density in the undisturbed stream. 

For the present problem of a wing performing small harmonic tor- 
sional oscillations of amplitude uq about some spanwise axis xq 
and small harmonic vertical translations of amplitude hQ (see 
fig. 1(b)), the equation for is 



where os is the frequency of oscillation. Substitution of eqiiatlon (4) 
into equation (2) gives 

w(x,y,t) = h + Ua + ^x - x^ja (5) 

Since equation (l) is linear, the velocity potential satisfying equa- 
tions (1) and (5) may be regarded as the sum of three potentials, each 
being associated with one of the terms on the right-hand side of equa- 
tion (5)» Thus the potential is obtained in the form 

( 6 ) 

Thickness effects are considered negligible and, as a result, 

the velocity potentials in equation (6) are associated only with lift 

and are antisymmetrical with respect to the plane Z = 0. Thus only 
one surface of the projected wing need be considered. The top svirface 
(Z = +0) is chosen in this analysis and, since the antisymmetrical 
potential is simply opposite in sign on the bottom surface, the pres- 
B^ire difference supported by the wing is obtained by means of equa- 
tion (3) as - . 
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EQUATIONS FOR VELOCITY’ POTENTIAL 
Integral Form 


The boundary-value problem for the arrowhead wing with supersonic 
edges, discussed previously, is similar to that for the wing of infinite 
span treated in reference 1. The problem in reference 1 is shown to be 
satisfied by a distribution or superposition of sources over the upper 
wing surface. Thus the velocity potential for the arrowhead wing, 
irrespective of the type of time variation, is 

ai (8) 

where w( |,ti) represents the surface variation of source strength and 


M(x - I) T R 

Tn ^ = X 


' 1,2 




R = ^(x - - p^(y - Ti)‘ 


The region of integration S, in equation (8), is the part of the wing 
(in I T] -plane) cut out by the reflected Mach cone opening upstream of 
the point (x, y,+0) as shown by the shaded area in figure 2. 

For the present problem of the wing performing small hapmonic 
torsional and vertical oscillations, the time variation of source 
strength is given by 

lent 

w(t) = e 

Thus, as shown in reference 1, equation (8) takes the form 


9((x,y,+0,t) 



~wU>Ti) -i(x-|)co 
R 


cos 


H) 


dT] d5 


where 


(uM _ 

2 ~ 2 

up 


( 9 ) 
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and the values of— ■ w(I,ti) for the harmonic case from equation (5) are: 
For ii, 

w( I ) = ItDho = — ho 

For Ua, 

w(§,T]) = UttQ 

For - Xq) , 

w(^,Ti) = (1 - XQ)ttQ 

For certain parts of the region of the wing (shaded in fig. 2) to 
which equation (9) applies, the Integration -cah be carried out as in 
the case of the infinite wing considered in reference 1. However, the 
entire potential 0 does not appear to be obtainable in terms of known 
functions. Previous mention of this point was made in the introduction. 
In order to give consistent treatment to the entire region of integra- 
tion, the integrand of eqijation (9) maybe expanded in Maclaurin's 
series with respect to < 0 , with the result that each term of the expan- 
sion can be readily integrated. Thus equation ( 9 ) is used in the 


following form: 



( 10 ) 

where [m/0 denotes the Integral part of m/2. 


Integrated Form 

The velocity potential to the third power of o) is considered 
sufficient for a large number of practical applications. Thus equa- 
tion ( 10 ) to the third power of 55, with any 55 that may come 
from w(5,ri,t) temporarily neglected, is 

0(x,y,+O,t) = - 

a^ -^ + bQR + b^^RjdTi d| 


v(|,Ti)(ao I + ai I + ag ^ + 


( 11 ) 
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W 


where 


_ o5r p 

= 1 - iCDX - + 1 


5^ 


x3 


= to + 


—2 

oa^x 


a>3 

i — 
^ 2 


ag 




X 


a 


3 


-1 


ns^ 

T' 


h 


0 


5^ 53^ 

2m2 ^ ^ 2m2 "" 



The question may he raised at this point as to whether equation (ll) 
still represents a solution of the houndary-value problem under discus- 
sion. With regard to this question, it can be shown that, when all the 
terms involving o) up to a given power are taken into account, the dif- 
ferential equation (l) is satisfied to this power. The boundary condi- 
tion of tangential flow, equation (5), is satisfied exactly, regardless 
of the order of o3 considered. 


Putting the values of w(^,ti) given in equation (9) into equa- 
tion ( 11 ) and performing the indicated surface integration (over the 
shaded region in fig. 2) yields the following forms for the velocity 
potentials in equation (6): 


K = - 


. 2 „2 2 „ -1 X + . T , -1 X - p^ly 

A^x - u y + Bi CCS cos . ') 


= - 1 




> ( 12 ) 


= 


. . da cos-^ . D. ccs-^ - 

XoFi(x,^ 
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The coefficlent-s A, B-[_, B2, C, and D2 are given in appendix A. 

In order to perform the integrations that led to the results represented 
hy equation (l 2 ), the use of characteristic coordinates (Mach lines in 
fig. 2) was found convenient. Reasons for using these coordinates are 
given in appendix A, where a brief account^f the derivation of equa- 
tion (12) is presented. 

In the limiting case where the Mach lines coincide with the leading 
edges of the wing, that is, when PX = 1, equations (l2) tsie the fol- 
lowing form: 

'^h = - i 

= - ^ Ai\/x2 - p2y2 

y 

The coefficients A]l *^1 ^hd the derivation of equation (13) are 

also given in appendix A. 


r (13) 


FORCES AND MOMENTS 


The preceding result” for the velocity potential, equation (l 2 ) , is 
now used to obtain the forces and moments on any one of a series of 
wings generally referred to herein as arrowhead type. As stated in 
the introduction, the wings being considered are those that may be 
formed from the delta or triangular wing by cutting the rearward part 
of this wing in such a manner that the resulting trailing edges lie 
ahead of the Mach cones emanating from their foremost points. Sketches 
of several wing plan forms having the aforementioned characteristics are 
shown in figure 3. 

The forces and moments acting on a wing section, such as section y 
of- figure 1 , are derived in general form, that is, in the form applicable 
to any of the aforementioned wings. These forces and moments are useful, 
for example, in a strip flutter analysis. Total forces and moments are 
also derived, but only for the delta wing. The damping part of the 
total moment due to aQ is used in a study of the possible loss of 
aerodynamic damping (loss Indicates that single-degree torsional flutter 
is possible) of an oscillating delta wing, which infers the same pos- 
sibility with regard to the other wings considered. 
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In deriving the expressions for the section and total force and 
moment coefficients, it is convenient to employ the variables p, y, 
and [iQ as the nondimensional quantities, obtained by dividing the 
variables x, y, and xq by the chord 2b of the midspan section of 

bcD 

the wing, and to introduce the reduced frequency It = — . 


Section Forces and Moments 


The local pressure difference between the upper and lower surfaces 
of the wing- is given by means of equation (?) as 

Ap = -2p 

The expression for the section force, positive downward, is therefore 


U_ 

2b 


(14) 


P = 


-2b 



Ap dp 


(15) 


and the section moment, positive nose up, about 
rotation x = Xq Is 



the arbitrary axis of 


dp 


(16) 


The limit of Integration p^ in equations (15) and (l6) has the fol- 
io-wing values with respect to the different plan forms shown in figure 3: 


For plan form A, 


For plan form B, 


For plan form C, 


For plan form D, 




^^1 


= 1 


“i = 1 ; 


i 




for 0<v 

m^_ m]_ + m2\2b ' 




I y_ 

2b m2 


for _ i] < V < 

m^ + tuq y2b y = — 


s 

2b 


(IT) 
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After equation (12) is substituted into equations (l5) and (l6) and 
tbe indicated integrations are performed, the results can be .written as 


, 2 2 icot 

P = -4pbU k e 




~(Li + cto(^3 + 


(18) 


and 


= -kpb^k^e^"*^ 


+ iMg) + cco(M3 + iMj 


(19) 


In equations (l8) and ( 19) the reduced frequency k, on which the L^’s. 
and Mj_ '_s are dependent^ is related to oi and cB in the following 
manner: 


k 


bto 
U ^ 




A brief account of the derivation of and expressions for the force and 
moment-coefficients and Mj_ (where i = 1,2, 3, 4) appearing in 

equations (l8) and (19; is given in appendix B. 

In equations (l8) and (19) the coefficients + iLg , + il<^ 

and L3 + iLi,. , M3 + are the section lift and moment coefficients 

associate'd with vertical translation and rotation, respectively. The 
real-component-of L;j^ + iL2 , for example, is in phase with the dis- 
placement ho and the imaginary part is 90° out of phase with this 
displacement. Similar interpretations can be given to the remaining 
component's. The imaginary part of each of the coefficients is propor- 
tional to the aerodynamic damping force or moment associated with the 
respective motion. 


Total Forces and Moments for Uelta Wing 


The expression for the total force, positive downward, on plan 
form A of figure 3 is 


P 



v/x 


Ap dp dv 


( 20 ) 
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Tlae total moment, positive nose up, about the arbitrary axis of rota- 
tion X = Xq is 


H 


a 



XQ)dM- dv 


( 21 ) 


After equation (12) is substituted into equations (2Q) and (21) 
and the indicated integrations are performed., the results can be written as 
(using barred letters throughout to designate total forces and moments) 


( 22 ) 


P = -8pbWe^^^ 


b 


+ H12) + cCq^L^ + H4) 


and 


= -Spb^uSc^e^"*^ 






(23) 


The derivation of the total force and moment coefficients and Mj_ 

(1 = 1,2,3, 4) is also given in appendix B. 


SOME CALCULATIONS AND DISCUSSION 


It may be of interest to the reader to examine the spanwise dis- 
tributions of the various components of the section lift and moment 
coefficients for a particular case. Thus, the components of equa- 
tions (18) and (19) ^ given more fully in equation (Bl), have been 
evaluated at different spanwise positions y for plan form A of fig- 
ure 3 for the following set of conditions: X = Hq = 0.5, 

= 1.75; and k = 0.04. These sample results are plotted as functions 
of the span position in figure i)-. The spanwise variations of the 
section-force components are shown in figure U(a) and the corresponding 
variations of the section-moment coefficients in figure ^(b) . 

The components of the total force and moment coefficients (given by 
equation (b4)) have also been evaluated for the previous set of condi- 
tions. The results, after being referred to the nondimens ional span of 
the wing 2X, are represented by dashed lines in figure 4. As should 
be the case, the areas under the dashed lines are equal to those imder 
the respective section-component curves. 

As a result of the linear dependence of the total force and moment 
coefficients on the, semiapex angle of the wing, as may be noted in 
equation (b4), the average ordinates of the distributions given in 
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figure 4 -will not change If- the apex angle ia changed. Only the shapes 
of the section force, or moment distributions ■will be affec-ted by such a 
■variation in angle. 

In_the example plotted, the total component of moment coeffi- 
cient M4 is negative. This term ■would therefore not contribute to 
the aerodynamic damping but would act as a source of energy for the 
oscillating system. This result is significant since it indicates the 
possibility of single-degree torsional flutter. 

The ■wing plan forms discussed in this paper exhibit the possibility 
of loss of aerodynamic damping in torsion for certain ranges of Mach num- 
ber M and axis-of=rotation location pq* This possibility is indi- 
cated, as mentioned in the previous paragraph, by the negative ■value of 
the torslonaL damping-moment component M4 in the example given. Since 
a wing oscillates as a rigid body at a fairly low frequency of oscilla- 
tion, the main results of- the loss-of-damping phenomenon can be obtained 
by maintaining H4 to the order of l/k only. Thus, there is obtained 
for plan form A of figure 3 the following result from equation (b4) : 

M4 = -^[T2 (m^ - l)p ^ - 4p (4 m^ - 5) + 3(SM^ - ^ (24) 

_In general, torsional stability depends on the sign and magnitude 
of Positive ■values of-'M^ indicate stability and negative ■values 

Indicate possible instability. The borderline case is thus given 
by M4 = 0. 

The ra^es of- values of- Mach number M and axis of rotation Pq, 
for which M4 in equation (24) vanishes, are shown in figure 5» Only 
one curve is obtained for all delta wings with supersonic edges, that 
is, PX ^ 1, beca\ise M4 is dependent on X for magnitude only, as 
may be noted in equation (24) . The region inside the curve in figure 5 
is the region of possible instability. Por convenience, a second 
ordinate is given in figure 5 showing the values of — X below which the 
leading edges become subsonic. With the aid of this second ordinate it 
may be seen that a delta wing with a semi vertex angle of 45° is the 
narrowest triangular wing that may show a torsional instability and yet 
retain the characteristic PX > 1 because the uppermost part of the 
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curve in the figure corresponds to a value of 1, = 1 when pX = 1. 

For a narrower wing to have at least sonic edges ^ it must he moving at 
a Mach number which would place it in the stable region of figure 5. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., July 11, 1951 
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APPENDIX A 

COEFFICIENTS OF VELOCITY-POTENTIAL EQUATIONS 


A derivation of equations (l2) and (13)__1 b given briefly here. 
The surface integral represented by equation (ll) ■when ■written in 
terms of the characteristic coordinates u,v becomes 



w(u,v)aj^ 


(u + v)^ du dv 

- ^)(^o - 


+ 



■w(u,v)bja(u + 


UQ - 11 




(Al) 


•where 

u * ||( I - 0T) ) 

Uq . ^ (x - py) 

V = (S + Pn) 

The limits of integration in equation (Al) , from the characteristic 
coordinates previously defined and from figure 2, are: 


Region 

^1 

Ug 


Vg 

I 

0 

^0 

0 

"^0 

II . 

1 

— V 

7 

^0 

7Uo 

0 

Ill 

yv 

0 

0 

^0 


(A2) 




3 


NACA TN 249.ll- 


17 


where 


_ 1 - PX 

' ~ 1 + px 

Putting the appropriate values of w(u,v) and the limits of integration 
given in equation (A2) into equation (A1) and performing the indicated 
integrations, and then reverting to xy-coordinates, yields the potentials 
in the forms given in equation (l2). Retaining only those terms 
associated with h or a that Involve <x> up to the third power yields 
the following values for the coefficients in equation (l2): 


6M2a2C 


M^(4c + 3) + a 






+ ( 3M^ + a) 




72 M^c 


M^(l8a^ + 31ff + 15) + 3<j(5cr + 3^ - 


1 1^ M^(38a -h 45 ) + 3ct(4<t -f 9 TI 

72 JTcr^ L 


2 2_2 

B,(x,y) . ^ (XX . y) - 1 ^ (XX . y)= - - 3) 

M^(2a + 5) + 3^ (Xx + y)^ 


(A3) 


(Xx + y)^ + i - 7 /p 

48M^a7/2 


B2(x,y) = B]_(x,-y) 


, M 3) . 1 + 19a h- 15) * 

^ 3a^ 24M^3 1_ \ 


cr(a + 3) 


Xp^xy^^ 


241^0^ 


M^(26c + 45) + 90 


■ y 

Continued on next page 
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Di(x,y) = - 


(Xx + y) - + y) + -^|^(Xx + y)^ - 




2ct' 


i + y)^ + 1 + 3)_(Xx + y)^ + 

^ ^/p + 5) - aCa - 3)1 (Xx + y)^ + 

km a' ^ -■ 

+ 5) + 3<7 (Xx + y)^ - 


X3^as^ 


(2a + 5)(Xx + y) 


%(x,y) = D^(x,*-y) 


J 


(A3) 


where 


2 2 

a p^X^ - 1 

For the case in which the velocity coln|jofients normal to the 
leading edges are sonic or in which the Mach lines coincide with these 
edges (that is, when pX = 1, 7 =s 0> or d c O), it may he noted from 

equation (A2) and figure 2 that-regions II and III no longer exist. 
Therefore, the potential in this case is determined by integration 
over region I only. The coefficients of equation (I 3 ) are thus found 
to be; 


. 2 1 . x35 

“ pf- ‘ ^ 3 *■ 9 om21: 


5 ^ P^y^ ) + (tx^ + 2p2y2j 




xoy 


210M^L 


7(x^ - pV) + M^(3x^ + 2 pV) I 


Ha4) 




315M- 


The use of characteristic coordinates leads to very simple expressions 
for the limits of integration given in equation (A2) . In addition, 
their use simplifies the derivation of equation (Ak ) . 
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APPENDIX B 

FOPCE AND JiKDMENT COEFFICIENTS 
Section Force and Moment Coefficients 


The coefficients Lj_ and {-where i = 1,2,3?^) are obtained 

hy substituting equation (l2) into equations (l5) and (l6) and grouping 
the results in the form given in equations (l8) and (l9) • These 
coefficients are therefore defined as follows: 


Li + iLg = L]_' 

+ iDi|_ — * 

+ iMg = M|_’ 

M3 + iM^ = M3 ’ 

•where 


+ iLg' 

+ iLi^'- + 2po)(Li« + iLg’) 

+ iM^* - + IL2*) 

+ 1M4* - 2po(L 3* + ILi^') - + 2^o)(Ml + IMg) 


Li’ + iLg' 



F]_(2u,2v)du - i F]^(2|i]_;2v) 


(Bl) 


L3' + ILi^* 



Fg(2u,2v)du - i Fg(2ui,2v) 


M]^’ + IMg' 



F3_(2u,2v)u lU - 2 U]_Ft_(2p3_,2v) 

_ 


Fg(2u,2v)u f 2U3_Fg(2u3_,2v) 

In equations (Bl) the quantity Fx(2u,2-v) is obtained from the qxxan- 
tity Fp(x,y), defined in connection with equation (l2), by replacing x 


M3' -KiMi,' 


4 

L^'vA 
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by 2|i, y by 2v, and regarding tbe 55 associated with the coeffi- 
cients kf and B2 in F]_(x,y) as 55' where 


55* = 



The quantity F2(2 |j., 2V) is obtained in a similar manner from the 
expression 


= C \/x2 - pV-+ co=- 1 + Ds coa-1 

the right-hand side of which is recognized as part of in 

equation (l2). “ 


Consider now the following general definitions, particular forms 
of which are contained in and (i =-1,2, 3, 4); 




e e+m+n 


F =2 
mn 


/«2. iH/, \n -1 1^1 + 

(s xv) m + v) COB ^ + 


(-^XV)« - V)-' COB- 


G =2 
mn 


'v/\ 


H 


,m+n+2 ^ r^l .n ./.2 „2..2 


H \/li - p V du 


r (B2) 


g e _ gS+m+n+l 
mn 




'v/x 


Ax. . v)“ OOB- d. 


,27-. 


J 


With the aid of-equation (B2), the L^' and M^' (i = 1,2,3, 4) of 
equation (Bl) may be written as 
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-2irL]_ ' 

-21CL2’ 

-2.L3' 

-2rfLi^' 

-2rtMj_' 

-2«M2’ 

-2rtM3' 

-2jtMi^’ 


■*■ ^2^02°] ^“2^02 ■'■ “3^20 ■’■ “4®03 


“5 " ®21 1^3^03 


° + ^4^04°)' 


■ J + (“1°01 ■ “2®02 ■ “3®20 ■'■ ^2®02° " P3®'o3°)’^ 


^1*^01 ^2®02 ■'■ ”^3^20 Vll ■’■ ^2^01 ' 


0 0 0 
63H02 + 54F2^2 + ^5^03 


■(^1®01 ^1^11 ®2^0l' ■*■ ^3^02 


”)^ • { 


^2*^02 ^3^20 


° + S H ° 

4 03 '5"^21 "4"'12 5 03 


+ 8,H_ + - 8gF._" - 6_F_ ' - 8„F.,''llsl 


*^-6F '-8F ^ 

13 7 03 8 04 


)' 


“2^03 ■'■ 
+ 


|“l(®02 ■ ^01] + Pl%l + ^2(^12 " ^2 ), 

°^3*^2l “4(^04 “ ^03) °^(^22 ~ '^21) ■'■ ^3%3 

P4(fi 4° - Hoi,°^ 

Pl(®bl " ^11°) k ■'■ [^1^02 ■'■ *^(^02 ■ ^03) ■'■ “3 (*^20 " ^21) 

P2Hi2° + P3(Hq 3° - Fi3°)k 

^1^02 ■*■ ^2 ^^03 " '^02) ■*■ ^3^®21 “ ^20) ■*■ ®1^21*^ ^2^11' 

®3Hi2 + ^ 4(^22° " %2°) + ^5(^13° ■ ^03°) 

n^Goi - E02 ) + - F2i°] + 52^Hoi’ - F^;l') + 

^3(^02 “ ^12 ) ^ [72^03 ’'3^21 ■'■ 71 ).^Go 3 " ® 0 l^ ■*■ 

75^021 - ^22 + ^k^22 + ^5%3° + ^6(^13° - F23°) + 
^7(®b3' - ^13 ’) + ^8(Ho4° - Is: 


HB3) 


_y 
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where 




“2 


“3 ■ ^ 2 ( 3 “® + “) 


+ Cf 




72 P ■- 


:i^(l8a^ + 310 + 15) + 3(j(5a + 3]] 


ttc 


M 


720 V 


^[^(38(1 + -Its) + 3o(li^+ 9)] 


P2 » --42 


3 ^ ^ 

^ 120 ^ 


^ [M 2 ( 2 a + 3 ) + o_ 


P4 


m\ 


480 ^a 


27 f 72 


M^( 2 a + 5 ) + 3 ^ 


’'l = a 


’'2 


M^X 


M^X 

73 = - 2 

■J a » 
M^X 


30 


+ 3) 7k = ■*■ ^5) + ^ 3)j 

M^(26a + 45) + 9"? 


M^X 

^^5 = "2402ct3 
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The coefficients in equation (B3) ^ that is> y-^, and 

originate > respectively^ from the coefficients A, B]_ or B 2 , C, 

and D]_ o't' D 2 defined in appendix A* The quantities and 

contained in equation (B 3 ) are easily determined from equation (B2) . 
The q-uantltles and however^ are detennlned hy integration^ 

as indicated in equation (B2) . The integral results that are required 
in equation (B3) arei 

_ 8/2 .2 2\3/2 

^01 = 3^^^1 ‘ P V j 

Gq 2 = 2 - p^v^) ^ - p^^ cosh"^ ^ 

O03 - 

G20 = ‘ coah*’- h 
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„ 0 1 „ 0 2 , -1 
HqI = 2X ^02 - ^r- ^ PV 


“02° = ^ “ 03 ° - ^ 


3X 




“03° = E “04°. - ^ 


- P^v^ + v^CScr + 5) cosh' 


.1^ 

3v 


TT ,0 ^ J:. T 0 2^^ 

5X ^05 15X 


1/2 


TT ® ^ •pi ^ J . -*~ -p 

^11 - -g^2 ^03 + 21 ^12 3X 


X.v^(2cr + 5) + X,3 jj,^2^2 ^^x^ 2 _ p 


2 2 
V 


1 „ 0 80^/2 ^2 ,/~2 „ 2.,2 


H' 


0 . _L TT 0 ill 

1 + r T Q • 


1/2 


12' = “7772 ^04 ^ ^3 - 35^2 


12X‘= 

v^(5cf + 3)cosh"^ 


5i2v2^x^ ^|iI^2T7V + 


1/2 r 


V V - ^ (220 4. 25) IV^ -H 
20 a» LyK ‘ 




H . 


a. 

12x3 

^04 ■ 

3i2 

^13 

v\0 

+ 3)cosh ^ 

PvJ 

1 

0 

1 

( 

^14 

30x3 

^05 “ 

6x^ 

Xv ^(60 + 5) 



1^0^/2 


3i2 1_ 


,2 2 rr" 

1 V 


2 2 

p V + 


A/2r 


3^ ^3 - ^5,3 


3 2 2 


31 V |J.^ + 
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% 




1 3^ „ 0 23^ „ 0 ^ „ 0 

^11 = "lx + -y ^13 " 2 ^22 - 

+ V^cosh"^ 


H 03 I - ^ - i^sV^[T3v4i2 . xv3(ac . 1 ) 






As may “be suspected from tlae complicated results generically represented 
by and 'B^, the section lift and moment coefficients 

are not simplified to any large extent by algebraic combination. The 
lift and moment coefficients, as given in equatipn (B3), are in a form 
which is perhaps better for numerical evaluation than a combined form. 

Some cautionary remarhs are in order with regard to equations 
and (B3) . In evaluating the coefficients Lj_ and Mi (where i = 
the teims involving should be ignored because the lift and moment 

coefficients associated with a have been treated completely only to 
the same order as those associated with h/U. In equation (B3) only 
real quantities are admitted. Therefore, the following restrictions 

must be noted (for ~ > iV 

\ ^^l P/ 


\jp^ 3^v^ = 0 

3(W+ v7 = ° 


Px - 

3(\ur-'^ ' ^ 


cosh-l y = 0 
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Total Force and Moment Coefficients 


As indicated in equations (20) and (21), the total lift and moment 
coefficients Li and Mi for the delta wing may be obtained from the 
previously derived Li and Mi by setting = 1 in Li and Mi and 

integrating the resulting expressions over the span of the wing. The 
results of the spanwise Integration combine_readily to yield expressions 
which are much sin^ler than the corresponding ones in the section- 
coefficient case. Without further detail, the total force and moment 
coefficients are: 


L-j^ + 1L2 — L^^ + ILg 


L3 + iLj^ -13+ iLj^ 


Mi + iMg * Mi + iM2* 


i V=- • - « 


- SpQ^Li + iL2 j 




M3 + iMi^. = M3' + IM^' - 2 hq(l3 ' + ILi^'j - (I + + iM^ 


J 


(b4) 


where 


-'2 = 


"3 

^4' 


M-, ' 


2X 

M^X 


' 153 '^ 

\ 


3k ■ 


2X 


33^ 


4x 

2M^X 

33k 

53^ 


(4m^ + 1)1 


■)' 
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4X , 

3P^‘5p5 ^ 





l6x 

1513^ 


Pk 



k 


If equations (22) and (23) had been put in the conventional form, 

Xj ~ C^qS 


M = CjjqSc 

where S Is the area of the wing (for the triangle S = 4b\) , q is 
the dynamic pressure, and c is the wing chord, then the lift and 
moment coefficients Ql and Cm, obtained from equation (b 4) , would 
he independent of X. This result is rather remarkable since it means 
that the lift and moment coefficients for a triangular wing with 
supersonic edges are independent of the vertex angle of the triangle. 
The same observation was previously noted in reference 4, 
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(b) Section y. 

oordlnate system and the two degrees of freedom a and h. 
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Figure 3*- JJifferent plan rorma ror whicli the force equation (18) and the 

moment equation ( 19 ) apply. 
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FigTire 5-- Curve showing ranges of M 6ind X - — against for which 

D U 

H 

the aerodynamic torsional damping moment vanishes for wide delta wings 






